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Abstract 

We develop a complex scaling method for describing the resonances of deformed nuclei and 
present a theoretical formalism for the bound and resonant states on the same footing. With 
31 Ne as an illustrated example, we have demonstrated the utility and applicability of the extended 
method and have calculated the energies and widths of low-lying neutron resonances in 31 Ne. The 
bound and resonant levels in the deformed potential are in full agreement with those from the 
multichannel scattering approach. The width of the two lowest-lying resonant states shows a novel 
evolution with deformation and supports an explanation of the deformed halo for 31 Ne. 
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I. INTRODUCTION 



The investigation of continuum and resonant states is an important subject in quantum 
physics. In recent years, there has been an increasing interest in the exploration of nuclear 
single particle states in the continuum. The construction of the radioactive ion beam facilities 
makes it possible to study exotic nuclei with unusual N/Z ratios. In these nuclei, the Fermi 
surface is usually close to the particle continuum, thus the contribution of the continuum 
and/or resonances is essential for exotic nuclear phenomena lM5|. 

Several techniques have been developed to study the resonant states in the continuum. 
One of them is the R-matrix theory in which resonance parameters (i.e. energy and width) 
can be reasonably^determined from fitting the available experimental data [gj]. The extended 
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have been also developed. The conventional 



scattering theory is also an efficient tool for studying resonances. More precisely, the scatter- 
ing phase shift method is adopted and the resonant state is determined from the pole of the 
S matrix qJ. Computationally, it is desired to deduce the properties of unbound states from 
the eigenvalues and eigenf unctions of the Hamiltonian for bound states so that the methods 
developed for bound states can still be used. For this purpose, bound-state-type methods 
have been developed, including the real stabilization method (RSM) 10j, the complex scaling 
method (CSM) 11|, and the analytic continuation in the coupling constant (ACCC) method 



121 ] . These bound-state-like methods have gained a great development. For examples, many 



efforts have 
RSM 
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Deen made in order to calculate more efficiently resonance parameters with the 
151 ] . Its extension to the relativistic framework has been presentes in Ref. 16| . 
where the resonance parameters obtained are comparable with those from the relativistic 
ACCC calculations. Combined with the cluster model, the ACCC approach has been used 

. An at- 
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to calculate the energies and widths of resonant states in some light nuclei 
tempt to extend the ACCC approach to the relativistic mean field (RMF) model was made 
in Ref. 19|. Further development of the formalism was presented in Ref. j^oj], where the 
wave functions of the resonant states were determined by the ACCC method. Based on the 
relativistic extension of the ACCC method, the structure of resonant levels was investigated 
and good pseudospin symmetry was disclosed in realistic nuclei 



221 ] . Although it involves 
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the solution of a complex eigenvalue problem which causes some difficulties in practice, the 
CSM has been widely and successfully used to study resonances in atomic and molecular 
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systems 







]23fl and atomic nuclei 
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was first outlined by Weder 



27] . Its generalization to the relativistic problem 



29j has 



Recent progress includes the following: Alhaidari 
presented a general and systematic development of an algebraic extension of the CSM to 
;he relativistic problem by expanding the Dirac spinors in Laguerre functions , Bylicki et al. 



30] have formulated a positive-energy-space-projected relativistic Hylleraas configuration- 



interaction method based on the complex coordinate rotation , and we have developed the 
CSM within the RMF framework in satisfactory agreements with the RSM, the scattering 
phase-shift method and the ACCC approach 31]. 

All these studies are mainly for spherical nuclei. Recently, the resonances of deformed 
nuclei have attracted additional attention. The interplay between the deformation and the 
low-lying resonance is very interesting for the open shell nucleus close to the drip line. The 
recent experimental data on the Coulomb breakup of the nucleus 31 Ne 32] have been well 
interpreted in terms of deformation with particle populating in the resonant levels The 
deformed halos have been investigated systemically in Ref. 5], where the resonances in 
the continuum play an important role. Some techniques have been developed to investigate 
the resonances of a deformed system. In Refs. 33|, |34|, the resonances as a function of 
deformation in an axially deformed Woods-Saxon potential and a quadrupole-deformed fi- 
nite square-well potential without a spin-orbit component have been studied by solving the 



coupled-channel Schrodinger equation. In Refs. |35|, |36], the ACCC method and RSM have 
been extended to investigate the resonances in axially deformed Woods-Saxon potentials. 
The single-particle resonant states in an axially deformed Gaussian potential without spin 
and tensor components have been investigated by the contour-deformation method in mo- 



mentum space 



371 ] . Using the multichannel scattering approach, Hamamoto 



38 



has studied 



how the single-particle energies change from bound to resonant levels when the depth of the 
potential is varied. This method has also been applied for the resonances of a Dirac particle 
in a deformed potential {39]. Hamamoto 40] has performed a systematic research on the 
resonances of deformed nuclei by the multichannel scattering approach, which includes the 
evolution of shell structure for exotic nuclei. However, as we know, CSM has not been used 
to investigate the resonances of deformed nuclei. Due to the great success of CSM in de- 
scribing the resonances for a spherical system, in this paper, we extend CSM to a deformed 
system, and examine its applicability and efficiency for the resonances of deformed nuclei. 
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II. FORMALISM 



Our purpose is to extend the CSM to describe the resonances of deformed nuclei. This 
scheme has a certain advantage, i.e., the bound and resonant states can be treated on 

'unctions for the resonant states are square 



231 ] . In the following, we start our scenario 



the same footing, because the complex scaled 
integrable just like those for the bound states 
with the single-particle Hamiltonian give as 

H = T + V, (1) 

where V is introduced with the axially symmetric quadruple-deformed potential consisting 
of the following three parts: 

Kent(r) = V f(r), 

Kou(r) = -f3 2 V k(r)Y 20 (#,cp), 

Ko(r) = -\vV Q g{r)(s-l). (2) 
Here k (r) = r^j^- and g(r) = — ^jp- - The parameter A is the reduced Compton wavelength 



of nucleon H/M r c. Similar to Ref. 38], a Woods-Saxon type potential is employed for / (r) 



\- . . To explore the resonant states, the CSM is used. A relative coordinate r in 

Hamiltonian H and wave function ip is complex scaled as 

U(9) : r -> re id . (3) 

Then, the transformed Hamiltonian and wave function are defined as Hg = U (9) HU (9) 1 
and ipg = U (9) if), where ipg is a square integrable function. The corresponding complex 
scaled equation is 

H ipo = E e tp e . (4) 



From ABC theorem 411]. the following is known: (i) a bound-state eigenvalue of H is also 
an eigenvalue of Hq, (ii) a resonance pole E — E r — iT/2 of the Green operator of H is an 
eigenvalue, Eg = E r — iT/2 of Hg; and (iii) the continuous part of the spectrum is rotated 
around the origin of the E— plane by an angle 29. 

To solve the complex scaled equation (j3J), it is convenient to adopt the basis expansion 
method. For the axially symmetrically deformed system, the parity it = (— 1)^ = (—1)' and 
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the projection of the total angular momentum along the symmetry axis Q are good quantum 
numbers; ipe can be expanded as 

^ = J>(0)&, (5) 

i 

where (pi = R n i(r)Yi mi ($,(p)xm. s (o'z) and the sum % runs over the quantum numbers n,Z,and 
m\ with Q = mi + m s . R n i(r) is the radial function of a spherical harmonic oscillator (HO) 
potential, 



*"<< r) = ^\j rwW L ' : - T (x2) e '"" " = 1 ' 2 ' 3 ' (6) 

x = r/b is the radius measured in units of the oscillator length b . Yi mi {&, ip) is the spherical 
harmonics. Xm s represents the spin wave function. The upper limit of the radial quantum 
number n is determined by the corresponding major shell quantum number N = 2 (n — 1)+/. 
Inserting the ansatz (J5J) into the complex scaled equation (j3J) and using the orthogonality of 
wave functions (pi one arrives at a symmetric matrix diagonalization problem, 

[T Vli + Vi^i] Ci = E Q c v , (7) 

i 

where T^i and V^j are presented as 

T, t = <r» UA-^-(l- + l±\ + J-) Ad f , (8) 



2M \dr 2 r dr/ 2Mr 2 
% = / fcV (re* 9 ) fadr . (9) 



Putting (pi into the above expressions, the matrix elements T^j are obtained as 

h 2 



Ti'i = e 



-i29 



2Mb 2 



^n(n + l + l/2)5 n /, n+ i + (2n + i - l/2)5 n> + ^(n - 1) (n + / - l/2)5 n> 

(10) 

Similarly, the matrix elements V^j are obtained as 

V£f = (<f>i>\ V ccnt (re ie ) |&) = V (n'l'\ f(re lf> ) \nl) 5 in S m , mi S m , ma , (11) 
Vff = (M Kou {re ld ) |&) = -f3 2 V (n'l'\ k (re ie ) \nl) (Z'mfl Y 2Q (0, ip) \l mi ) 5 m>ls , (12) 
V$° = (<f>i>\ V so \<f>i) = ~W (n'l'\ g(re ld ) \nl) (I'm^m'J (s ■ l) \lmim.) . (13) 

In these expressions, (I'm 1 ^ Y 2 q <p) \lm{) and (J.'m'pn 1 ^ (^s ■ tj \lmim s ) can be calculated by 
the usual angular momentum formulas, and the radial parts of integration can be calculated 



with the Gauss quadrature approximation. With the matrix elements 7V,t, V^f 1 *, and 
V^°, the solutions of the complex scaling equation (j3J) can be obtained by diagonalizing the 
matrix Hg. The eigenvalues of Hg representing bound states or resonant states do not change 
with 9, while the eigenvalues representing the continuous spectrum rotate with 9. The former 
are associated with resonance complex energies E — zT/2, where E is the resonance position 
and T is its width. 



III. NUMERICAL DETAILS AND RESULTS 

With the theoretical formalism, we explore the resonant states for the realistic physical 
system. In the numerical details, the complex-scaled Schrodinger equation is solved by 
expansion in the HO basis with 60 oscillator shells. The oscillator frequency of the HO basis 
is fixed as fkuo = 41t4 _1//3 MeV. For the comparison with the results from the multichannel 



scattering approach 



31 Ne is taken as an example. The corresponding parameters are 
4fl, i.e., a = 0.67 fm, V = -39 MeV, v = 32, and the 



chosen to be same as those in Ref. 
nuclear radius R varying with the mass number A of the system with R = r^A 1 ^, where 
ro = 1.27 fm is used. To achieve the information on the resonant states, we diagonalize Hg. 
The eigenvalues of Hg for these states with the quantum numbers varying from Vt v = ^ 
to = | are plotted in Fig.l, where the rotation angle 6 = 18° and the quadruple 
deformation 02 — 0.1 are adopted. All the eigenvalues of Hg, which correspond to the 
bound, the resonant, and the continuum, are respectively labeled as open squares, red open 
circles, and open circles. The solid line marks the position of the continuum with the rotation 
angle 29. From Fig.l, one sees clearly that the eigenvalues of Hg fall into three regions: the 
bound states populate on the negative energy axis, while the continuous spectrum of Hg 
rotates clockwise with the angle 29, and resonances in the lower half of the complex energy 
plane located in the sector bound by the new rotated cut line and the positive energy axis 
get exposed and become isolated. In realistic calculations, because the finite basis is used, 
the continuous spectrum of Hg consists of a string of points. 

Next, we demonstrate how the resonant states become exposed by complex rotation. The 
eigenvalues of Hg with 9 = 8°, 12°, 16°, and 20° are respectively displayed in Figs.2(a)-2(d), 
where the other parameters are the same as those in Fig.l. From Fig.2 one can see that, when 
6 = 8°, the resonances still have not been completely separated from the continuum, and to 
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identify the resonances in the continuum is still relatively difficult although several points are 
marked as the resonances in the complex energy plane. When the rotation angle is added to 
9 = 12°, the resonant states are clearly isolated from the continuum and distinguishing the 
resonances from the continuum becomes easy. When the rotation angle is added to 9 = 16°, 
the resonant states are completely separated from the continuum, and clearly exposed in 
the complex energy plane. This is because the continuous spectrum of Hg rotates clockwise 
with 9, while the resonances are almost independent of the rotation angle 9. From Figs. 2(c) 
and 2(d), the rotation angle is incresed from 9 = 16° to 9 = 20°, and the change of position 
of the resonant states in the complex energy plane is almost not observable. This indicates 
that the resonant states can be unambiguously determined at the same time as long as the 
rotation angle is chosen to be large enough, which is one of the advantages of CSM. In the 
present calculations, the exposure of the resonant states relies on the scale of the rotation 
angle, while their position in the complex energy plane remains almost unchanged with the 
variation of 9; i.e., the energies and widths are nearly 9 independent, which also reflects that 
the stability of the results with respect to the variation of 9 is considerably satisfactory. 

As we look forward to the resonances of the deformed nuclei states, it is interesting to 
observe the movement of the resonant states in the position with deformation in the complex 
energy plane, which is helpful to recognize intuitively the resonances of the deformed nuclei. 
In Fig. 3, we show the eigenvalues of Hg with deformation /3 2 = 0.0, —0.1, —0.2, and —0.3. 
The other parameters are the same as those in Fig.l. When the system is spherically 
symmetrical [fii = 0.0), only three resonances are observed, which correspond respectively 
to the resonant states I/7/2, 1/5/2; an d lgg/2- When the spherical symmetry of the system is 
broken (/3 2 7^ 0.0), the 3 resonances are broken down into 12 resonances, which correspond 
respectively to the resonant states l/2[330], 3/2[321], 5/2[312], 7/2[303], 1/2[310], 3/2 [301], 
5/2[303], l/2[440], 3/2[431], 5/2[422], 7/2[413], and 9/2[404] (denoted with the asymptotic 
quantum numbers Q[Nn z A\). The location of the 12 resonant states in the complex energy 
plane appears to be a significant movement when the deformation is changed. With the 
development of deformation towards the direction of larger deformation, the 12 resonant 
states in the position are more clearly separated. The similar movement of the position 
towards the prolate shapes has also been obtained, and not been displayed here. 

Although the resonances of several points of deformation are displayed in Fig. 3, the 
resonance parameters of these points have not been fully determined for the resonant states 
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exposed in Fig.l. In the following, we check how the results depend on 9 in the present 
calculations in order to present the accurate values of the resonance parameters. According 
to the mathematical theorem the resonance energies should be independent of 9 (for large 
enough 9). However, in the present numerical approximation, resonance parameters are 
dependent slightly on 9 due to basis expansion with finite shells, there will be no eigenvector 
that is completely independent of 9. The resonance energies will move along trajectories in 
the complex energy plane as a function of 9. The best estimate for a resonance energy is 
given by the parameter 9 value for which the rate of change with respect to 9 is minimal. 
To locate the minimal rate point, we perform repeated diagonalizations of the eigenvalue 
problem of Hg with different 9 values in every deformation concerned. The illustrated 
results are shown in Fig.4 for the resonant state 7/2[303] of the deformation point j3 2 = 0.3, 
where N represents the number of oscillator shells of the basis. From Fig.4 one sees clearly 
that, for different N, the dependence of the resonance parameters on 9 shows a different 
trajectory, and the range of the trajectory becomes smaller with the increasing N when 
N ^ 50. In every ^-trajectory, there exists a point corresponding to the minimal rate of 
change of the resonance parameters with respect to 9, which presents the optimal value of 
the resonance parameters. Despite the difference in the 9 trajectory, the optimal value of 
the resonance parameters tends to the same result with the increasing N, especially when 
iV ^ 50. Although the movement of the optimal value in the complex energy surface is 
observable, that is because we display these results in a small energy scale. Compared with 
the continuous spectra, the movement is negligible, especially when N ^ 50. Hence, the size 
of the basis being chosen as N = 60 is good enough in the previous calculations. To more 
accurately determine the value of the resonance parameters, a close up of the 9 trajectory 
is plotted in Fig. 5 for N = 60. In the vicinity of 9 = 17°, the resonance parameters are 
almost independent of 9, i.e., ^g- ~ 0, which presents the optimal value of the resonance 
parameters. 

By using the 9 trajectory, the resonance parameters of any deformation point can be 
determined for the resonant states exposed in Fig.l. In Fig. 6, we show the evolution of 
the single-particle energy with the deformation /3 2 for the bound and resonant states. In 
comparison with Ref.|4j, consistent results are obtained for the bound levels and the low- 
lying resonant levels. In addition, several higher-lying resonant levels are also obtained in 
the present calculations. Especially, the width is obtained for all the resonant states. In 
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Fig. 7, we plot the width for the four lower-lying resonant states 1/2 [330], 3/2[321], 5/2[312], 
and 7/2[303]. For the resonant state 7/2[303], the resonance energy and the width increase 
monotonously with the deformation varying from 2 = —0.4 to 0.6. For the resonant state 
5/2 [312], the energy and the width display an almost horizontal change with the deformation. 
For the resonant states 1/2 [330] and 3/2 [321], with the deformation varying from 2 = —0.4 
to 0.6, the energy decreases monotonously, while the width shows a unusual trend. Close to 
02 — 0.0, the value of the width appears to be at a minimum, which implies the spherical 
shape is more stable for the two resonant states. After the minimum appears, the resonance 
width increases with the increasing of deformation, which means these nuclei states become 
even more unstable as their levels become lower. With a continued increase in deformation, 
the width appears to be at a maximum value. Therefore, it is relatively difficult for the 
halo to forming in the deformation interval from 2 = 0.0 to the position of the width 
appearing at a maximum value for 31 Ne. Thereafter, the resonance width decreases with 
the increase in deformation. When the deformation increases beyond 02 = 0.25, the level 
l/2[330] becomes weakly bound. In the vicinity of 2 = 0.3, the levels l/2[330] and 3/2[202] 
appear to be crossing. With a further increase in deformation, the energy gap between 
l/2[330] and 3/2[202] gets larger and a maximum appears at 2 = 0.43. Thus, in this 
deformation interval from 2 = 0.25 to 2 = 0.43, the 21th neutron can occupy the weakly 
bound level 1/2 [330] or 3/2 [202], and the halo can be formed there. The most likely position 
for forming a halo should be near 2 = 0.3. This information indicates that 31 Ne may be a 
halo nucleus with deformation 02 ~ 0.3. 



IV. SUMMARY 



In summary, the complex scaling method is extended to describe the resonances of de- 
formed nuclei. A theoretical formalism is presented in which the bound and resonant states 
are treated on the same footing. 31 Ne is chosen as an example; the utility and applicability 
of the extended method are demonstrated. It is found that the exposure of resonant states 
relies on the scale of the rotation angle, while their position in the complex energy plane 
remains almost unchanged with the variation of 6, which reflects that the stability of the 
results with respect to the variation of 9 is considerably satisfactory. The movement of the 
resonant states in the position with deformation is exhibited in complex energy plane, where 
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it is seen that 3 degenerate resonances are separated into 12 resonances with the develop- 
ment of deformation towards the oblate or prolate shapes. The bound and resonant levels 
obtained for 31 Ne are in agreement with those from the coupling-channel calculations, where 
the resonances are determined by a multichannel scattering approach. With the deforma- 
tion evolution from the oblate to prolate shapes, the bound and resonant levels show a clear 
shell structure. Especially, the width of the two lowest-lying resonant states shows a novel 
evolution with deformation, and supports an explanation of the deformed halo for 31 Ne. 
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Fig 2: (Color online) The resonant and continuous spectra varying with the complex rotation 
angle associated with the model. Except for the complex rotation angle, the other parameters are 
the same as those in Fig.l. 
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Fig 3: (Color online) The movement of the resonant states in the position with deformation in 
the complex energy plane, where the development of deformation is towards the direction of oblate 
shapes and the parameters adopted are the same as those in Fig.l except for a different deformation 
in every configure. 
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Fig 4: (Color online) The 9 trajectories corresponding to the several different numbers of oscillator 
shells of the basis for the resonant state 7/2[303] of the deformation point fa = 0.3, where N is 
the quantum number of the main shell of the oscillator basis, and 9 varies from 8° to 22° by steps 
of 0.5°. 
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Fig 5: The same as Fig. 4, but a close up of the 9 trajectory in Fig. 4 for N = 60, where the arrow 
marks the position of the resonant state. 
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Fig 6: (Color online) Neutron single-particle levels in Woods-Saxon potentials as a function of the 
quadrupole deformation parameter /?2- The depth, the diffuseness, and the radius of the potential 
are -39 MeV, 0.67 fm, and 3.946 fm, respectively. Every level is labeled with the asymptotic 
quantum numbers f)[Nn 2 A]. 
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Fig 7: (Color online) The same as Fig. 6, but for the width of the four resonant states 1/2 [330], 
3/2[321], 5/2[312], and 7/2[303]. 
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